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O ; 1 (Pre-)History of the Model 



The Gonihedric 3D Ising model is a lattice spin model in which planar Peierls 
boundaries between + and — spins can be created at zero energy cost. In- 
stead of weighting the area of Peierls boundaries as the case for the usual 3D 
^ ■ Ising model with nearest neighbour interactions, the edges, or "bends" in an 

interface are weighted, a concept which is related to the intrinsic curvature of 
the boundaries in the continuum. 

The model is a generalised Ising model living on a cubic 3-D lattice with 
■ nearest neighbour, next to nearest-neighbour and plaquette interactions. The 

^ I ratio between the couplings of these three terms is fixed to a one parameter 

^ ■ family which endows the model with unusual properties both in and out of 

equilibrium. Of particular interest for the discussion here will be that the 
model manifests all the indications of glassy behaviour without any recourse 
to quenched disorder, whilst still possessing a crystalline low temperature 
phase in equilibrium. 

I In these notes wc follow a roughly chronological order by first reviewing 

Cf^ ■ the background to the formulation of the model, before moving on to the 

elucidation of the equilibrium phase diagram by various means and then to 
^\ • the investigation of the non-equilibrium, glassy behaviour of the model. We 

I apologize in advance for our narrow focus on things Gonihedric at the expense 

. of other lattice models with glassy behaviour, since the aim is to concentrate 

on giving an overview of the Gonihedric Ising model in 31?. 

The model has an unusual genesis since it was originally introduced as a po- 
I tential discretization of string theory. The Nambu-Goto [1] action (or Hamil- 

. tonian in statistical mechanical language) in bosonic string theory is given by 

' the area swept out by the string worldsheet as it moves through spacetime. 

Directly discretizing euclideanized versions of this action produced ensembles 
of surfaces, i.e. string worldsheets, which were dominated by collapsed and 
irregular configurations such as that in FiglT] and which were unsuitable for 
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Fig. 1. A typical collapsed triangulated surface resulting from a simulation of the 
Gaussian Hamiltonian in equ[T] 

taking any sort of continuum limit [2]. This was presumably a reflection of 
the well-known tachyonic instability of the bosonic string and therefore not 
entirely unexpected. The interest was in whether any straightforward modi- 
fications of such Hamiltonians did allow a continuum limit, and hence some 
insight into stable string theories in physical dimensions. 

Discretizing the Polyakov [3] form of the string action, which is equivalent 
to the Nambu-Goto formulation, by triangulating an embedded surface gives 
the simple Gaussian Hamiltonian 



where the vectors Xi live on the vertices of the triangulated surface and 
determine where the vertices sit in the embedding space. The sum {ij) is 
carried out over the edges of the triangulation and in principle one should 
also sum over different triangulations of the surface as a discrctized version 
of the sum over metrics in the continuum string theory with the Polyakov 
action. In any case, the observed behaviour of surfaces turns out to be rather 
similar in ensembles with and without this sum. One possibility for stabilising 
the surfaces produced by such Hamiltonians is to add an extrinsic curvature 
term as suggested first in the continuum by Polyakov and Kelinert [4J , which 
acts to prevent the surfaces collapsing 




(1) 
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H=- ^(X, - X,f + A ^ (1 - • n,) (2) 

where the rii^nj are normals on adjacent triangles Ai, Aj. The second term 
tends to align the normals on adjacent triangles and hence to "stiffen" or 
flatten out the surface _5j, at least over short length scales. The issue for Monte 




Fig. 2. A typical uncoUapsed triangulated surface resulting from a simulation of 
the Gaussian plus extrinsic curvature Ifamiltonian in equ[2]with A — 1.1. 

Carlo simulations is to determine whether this flattening is also effective on 
macroscopic length scales. 

Another possibility, which also uses a geometrical concept and which has 
some appealing features from the string theory point of view 6 , is to take an 
action proportional to the linear size of a surface - a notion first introduced 
by Steiner ^ . This can be discretized as 

(3) 

on triangulated surfaces, where 9{aij) = \tt — aij\ where is the dihedral 
angle between the embedded neighbouring triangles with a common link {ij). 
The d{aij) term also acts to flatten out the surfaces [5^. The elements making 
up these various possible surface Hamiltonians are shown in Fig[3] The origin 
of the word Gonihedric is a combination of the Greek words gonia for angle 
(referring to the dihedral angle) and hedra for base or face (referring to the 
adjacent triangles). 

Savvidy and Wegner posed the question of how to translate the essential 
features of such a Hamiltonian to a surface composed of plaquettes on a cubic 
(or hypercubic) lattice |9j. This was motivated by the observation that it was 
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Fig. 3. Two adjacent triangles in a triangulation of a surface with a common edge 
(ij), showing the normals n, the co-ordinates of the endpoints Xij and the dihedral 
angle aij 

possible to rewrite the Hamiltonian for an ensemble of such surfaces by using 
the geometrical spin cluster boundaries of a suitable generalised Ising model 
with nearest neighbour, next to nearest-neighbour and plaquette interactions 
to define the surfaces and their appropriate energies. It is generally much easier 
to simulate an Ising spin model than a gas of surfaces so one might expect 
some advantages, both numerical and analytical, from such a reformulation 
of the system. 

On a fixed cubic lattice all the edge lengths |X, — Xj\ will be identical, 
so the statistical weight of a surface configuration will be entirely determined 
by the 0{aij) = |7r — a^l factors, where the are restricted to multiples 
of tt/2 radians. In other words, the statistical weight of a plaquette surface 
configuration will be determined entirely by the number of bends and self- 
intersections it contains. There is no (bare) weight for the area of a plaquette 
in the surfaces in contrast to the usual 3D Ising model with only nearest 
neighbour interactions. 

The generalised Ising models which are used in the construction are of the 
form 

H = Ji'^aiaj + J2 ^ aiaj + J3 ^ CFiUjCrkcri. (4) 

{io) {{hi}) [i,j,k,l] 

where the Hamiltonian contains nearest neighbour next to nearest 

neighbour (((«, j))) and round a plaquette {[i,j,k,l]) terms. The couplings in 
such models can be related to the energy cost for a plaquettes in the Peierls 
interface a right-angled bend between two such adjacent plaquettes, Pci 
and the energy cost, for the intersection of four plaquettes having a link 
in common 

Pa = 2 Ji + 8J2 
Pc = 2 J3 - 2 J2 

Pi = -4 J2 - 4 J3 . (5) 

General cases of such Hamiltonians and their equivalent surface formula- 
tions had been studied in some detail by Cappi et. al. and found to have a very 
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rich phase structure for generic choices of the couphngs [10) , including homo- 
geneous, lamellar, disordered and bicontinuous phases. The Gonihedric model 
constitutes a particular one-parameter slice in this space of Hamiltonians: 

= -2k^o-,o-j + (Jtaj-^-^ arajCTkcri. (6) 

and for such a ratio of couplings /3yi = 0, which is the desired zero weight 
for plaquette area suggested by translating the Hamitonian of equ.®. The 
energy of the spin cluster boundaries for this Gonihedric model on a cubic 
lattice is simply given by 

E — n2 + 4^714, (7) 

where n2 is the number of links where two plaquettes on a spin cluster bound- 
ary meet at a right angle, ri4 is the number of links where four plaquettes meet 
at right angles, and k is a free parameter which determines the relative weight 
of a self-intersection of the surface. 

The particular ratio of couplings in equ.([6]) also introduces a novel semi- 
global symmetry into the Hamiltonian, as we discuss in the next section. This 
is related to a zero-temperature high degeneracy point where it is possible 
to flip non- intersecting planes of spins at zero energy cost. This symmetry is 
further enhanced at k = where the action becomes a purely plaquette term 

= -i ^ a.ajakcri. (8) 

and any plane of spins, even intersecting ones, may be flipped at zero energy 
cost at T = 0, which suggests a highly degenerate ground state for the k, — 
Gonihedric model. Series expansions and cluster variational calculations, 
described below, show that the degeneracy of the Hamiltonian is broken in 
the free energy for k 7^ 0, indicating a ferromagnetic low temperature phase in 
this case, but the nature of the order in the low temperature phase for /c = 
remains rather mysterious. 



2 Equilibrium Behaviour, by various means 

In what follows we use the standard definitions of quantities such as the 
magnetisation 

i 

and the Binder cumulant 

in discussing critical properties, where the () brackets denote thermal averages. 
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The energy {E) is just given by the average of the Hamihonian, the specific 
heat C by the variance of the energy and the susceptibihty x by the variance 
of the magnetisation. The scahng exponents are then defined by the singular 
properties of these quantities at the critical point (s) of the model 



where t = |(T — Tc)/Tc| is the reduced temperature, ^ is the correlation 
length and v is the correlation length exponent, giving the divergence of the 
correlation length at criticality. 

These can also be recast as finite-size scaling relations 



for lattices of linear extent L, which is more convenient when Tc is not known 
exactly, which is usually the case. 

Various methods have been used to investigate the phase diagram of the 
Gonihcdric model, ranging from generalised mean field techniques, through 
cluster variational calculations to Monte Carlo simulations. The gross phase 
structure of the Gonihedric model was apparent even in the earliest Monte 
Carlo work [TT], and is shown schematically in Figdl For n — Q and small 
values of k, there is a first order transition [12j . signified by a jump in the 
energy, E, at the transition point and first order (volume) scaling in the 
specific heat. 

As K was increased the transition softened to second order, but Monte 
Carlo simulations and other approaches have struggled to produce a consistent 
set of critical exponents. The first direct Monte-Carlo simulations [TT], for 
instance, found v = 1.2(1), 7 = 1.60(2),/? = 0.12(1) when k = 1, but later 
work using the scaling of the surface tension of a spin interface [13] obtained 
u = 0.44(2) and -f/u = 2.1(1). 



C = B + Cot-'^, M = Mot^, X = Xot-\ e = ^ot 



(11) 



C = B' + C^L"^, M = M[,L-, X = XoLi 



(12) 



K 



1 



2 order 




0.55 



Fig. 4. The phase transition line as k is varied. 
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Other approaches also gave results at variance with these: a low temper- 
ature expansion by Pietig and Wegner found a = 0.62(3),/? = 0.040(2) and 
7 = 1.7(2) at K = 1 pil; whilst Cirillo et.al. [H] found /3 = 0.062(3) and 
7 — 1.41(2) using a combination of the cluster variational method and Pade 
approximants ("CVPAM") [16]. 

These cluster variational calculations, ground state enumerations and 
mean field calculations all use a similar starting point - the decomposition 
of the lattice into elementary cubes and we now discuss these methods in 
some detail since they are useful in sketching out the phase structure. Indeed, 
the cluster variational method, allied with Pade approximant techniques has 
succeeded in providing a very plausible picture of the Gonihedric phase dia- 
gram by looking at an extended two-parameter Hamiltonian. 

One of the elementary cubes from the 31? cubic lattice is shown in FiglSl 
The full lattice Hamiltonian may be written as a sum over individual cube 
Hamiltonians 

= ^ X! "^^^J " i XI + X! Oia^a^ai (13) 

which differs from the full Hamiltonian by the symmetry factors in front of 
each term. If the lattice can be tiled by a cube configuration minimising the 
individual then the ground state energy density is — min h^. 
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Fig. 5. The ground-state and mean field analyses of the Gonihedric model require 
working with elementary cubes rather than single sites to capture the various pos- 
sible spin configurations. 

We list some of the inequivalent spin configurations on a single cube and 
their multiplicities in Table. 1. In the list of spins the first column represents 
one face of the cube and the second the other and the notation V'o,--- is 
borrowed from |10| The antiferromagnetic image of a configuration is obtained 
by flipping the three nearest neighbours and the spin at the other end of the 
cube diagonal from a given spin and is denoted by an overbar. 
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Table 1. Elementary cube ground state energies 



State- 


Top 


JJcrttom 


Eiicrgj' 


.Multiplicity 


i>o 


+ + 

+ + 


+ + 
+ + 


-3/2 - 3k/2 


2 


V'o 


- + 
+ - 


+ - 
- + 


-3/2 + 21k/2 


2 






+ + 
+ + 


-3/2 - 3k/2 


6 




+ + 


+ + 


-3/2 + 5k/2 


6 



With the Gonihedric vahics of the coupHngs there is a freedom to flip spin 
planes in the ground state as tpo, which would represent a ferromagnetic state 
when used to tile the lattice, and tpQ which would represent flipped spin layers, 
have the same energy for any value of k. The degeneracy increase when k = 0, 
the club of states of energy —3/2 is now composed of ■^o, V'O) V'6, V'g- 

This degeneracy means that the definition of an order parameter for the 
low-temperature phase is a moot point. The standard magnetisation M will 
not serve at zero temperature as is clear from consideration of the ground state 
energies above, since the freedom to flip spin planes means that it is identically 
zero. As we have noted, it appears that this symmetry is broken for k ^ 
at non-zero temperature, but in Monte Carlo simulations the magnetisation 
is still indistinguishable from zero, presumably because of the difficulty of 
removing such flipped planes within the timespan of the simulation. Even 
staggered magnetizations do not suffice as order parameters as the interlayer 
spacing between the flipped planes of spins can be arbitrary. 

One possibility on a finite lattice is to force the model into the ferro- 
magnetic phase with a suitable choice of boundary conditions such as fixed 
boundary conditions. Although this allows the use of a standard magnetic 
order parameter, it does some violence to both the scaling properties and the 
ground state structure structure itself. The corrections to first-order scaling, 
for example, are much stronger with fixed boundary conditions than with the 
more familiar periodic boundary conditions, although simulations of the k = 

Table 2. Scaling laws for Periodic versus Fixed Boundary Conditions for the critical 
temperature, maxima of the specific heat and susceptibility and minimum of the 
Binder cumulant. In all cases a d-dimensional cubic lattice of length L is considered. 





P.B.C. 


F.B.C. 


Cmax (-^) — 
Xmax{L) = 

Umin (-^) ~ 


/3e(00)+-^+0(^) 

70 +72i' + O(i^) 
$0 + ^+O{j^) 


/3e(oo) + ^+0(^) 

60 + CiL'' + OiL''-') 
Bo + ^ + Oi^) 
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model with fixed boundary conditions have given good agreement with this 
scaling [T7] . 

The approach taken to enumerating the ground states by splitting the 
lattice into elementary cubes can also be applied to mean field theory for the 
model and to extensions of mean field theory such as the cluster variational 
method. In the mean field approximation the spins are replaced by the average 
site magnetizations and an entropy term introduced to give the free energy. 
In the Gonihedric model the total mean field free energy may be written as 
the sum of elementary cube free energies </)(mc) 

4>{mc) = X/ "^'"^i + ^ X/ "^i'^j ^ ~~T~ ^ miUijmkmi 

+ ^ [(1 + mi}ln{l + m,) + (1 - m,)ln{l - m,)] (14) 
16 

iCC 

where mc is the set of the eight magnetizations of the elementary cube. This 
in turn gives a set of eight mean-field equations 

=0 (15) 

om, (j=i...8) 

(one for each corner of the cube) rather than the familiar single equation for 
the standard nearest neighbour Ising action. 

If we solve these equations iteratively we arrive at zeroes for a paramag- 
netic phase or various combinations of ±1 for the magnetised phases on the 
eight cube vertices, and the mean field ground state is then given by gluing 
together the elementary cubes consistently to tile the complete lattice, in the 
manner of the ground state discussion. Carrying out this program gives a 
rather simple mean field phase diagram for the Gonihedric model with action 
equ.®, with a single transition from a paramagnetic phase to a degenerate 
"layered" phase that is pushed down to /3 = at large k. The low tempera- 
ture phase is generically of the ■00,6 type, apart from k = where we see a 
tpQ g phase that is degenerate with these. A more sophisticated treatment us- 
ing the cluster variational method, described below, suggests that the layered 
phase has a slightly higher free energy than the ferromagnetic phase so the 
low temperature phase is in fact the ferromagnetic one when the model is in 
the regime with a continuous transition. 

The cluster variation method, or CVM for short, is based on a truncation 
of the cluster (cumulant) expansion of the free energy density functional on 
which the variational formulation of statistical mechanics is based 1181. Un- 
like mean field theory it generally locates rather accurately the boundaries 
between different phases in complex phase diagrams and, using the recently 
proposed cluster variation-Pade approximant method one can extract non- 
classical, precise estimates of the critical exponents. 

For the Gonihedric model, it is again possible to use the cube approxi- 
mation of the CVM, in which the maximal clusters are the elementary cubic 
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cells of our simple cubic lattice; for this approximation the free energy density 
functional has the form 



where Hg, is the contribution of a single cube to the Hamiltonian, C{x) — 
a; In a;, pa with a = 8 (4, 2, 1) denotes the cube (respectively plaquette, edge, 
site) density matrix, and the sums in the entropy part are over all plaquettes 
(edges, sites) of a single cube. 

Since the cluster variation method can be viewed as a generalised mean 
field theory, it is clear that it can give only classical predictions for the crit- 
ical exponents. In order to overcome this difficulty, one can use the cluster 
variation-Pade approximant method (CVPAM), which has proven to be a 
quite accurate technique for extracting exponents. The basic idea of the CV- 
PAM is that, since the CVM gives for Ising-like models very accurate results 
far enough from the critical point, one can try to extrapolate these results 
in order to study the critical behaviour. In order to determine the critical 
exponent of an order parameter Af, for example, one calculates M{(3) with 
the CVM up to a temperature at which the error can be estimated to be very 
small (typically of order 10~^), and then constructs, by a simple interpolation, 
Pade approximants for the logarithmic derivative of M{j3): the pole and the 
corresponding residue of each Pade approximant are then estimates for the 
critical temperature and for the critical exponent respectively [16j . 

The CVPAM calculations give a phase diagram whose topology agrees 
with Monte Carlo simulations [T3]. The value for the magnetic exponent is 
found to be /3 = 0.062 ± 0.003. Another important observation from the CV- 
PAM calculations is that at finite temperatures there is a violation of the flip 
symmetry of the Gonihedric Hamiltonian when k ^ 0, so that in the ordered 
region of the model the ferromagnetic phase is always stable with respect to 
the lamellar phase. 

Applying CVPAM methods to a slight extension to the Gonihedric Hamil- 
tonian gives useful insight into the nature of the phase diagram |19j . The 
Hamiltonian is modified to 



with a second parameter j which reduces to the Gonihedric model, with suit- 
ably rescaled couplings, when j = —0.25. The CVPAM calculations (and 
independent transfer matrix calculations [20]) for this modified Hamiltonian 
give a phase diagram whose k = 1 slice is shown in FiglBl 




(16) 




1 



K 



(17) 
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paramagnetic 




2 



lamellar 



ferromagnetic 



-0. 



1.25 



J 



Fig. 6. A K = 1 slice of the two-parameter Hamiltonian phase diagram. 

The structure of this phase diagram gives some strong hints as to why 
the critical exponent measurements in Monte Carlo simulations have been 
so inconsistent. The dotted line in Figgis the usual 3-D ferromagnetic Ising 
transition, along which the standard 3D Ising critical exponents such as = 
0.6294 would be observed. The line separating the lamellar and ferromagnetic 
phases appears to be bent slightly to the left according to both CVPAM and 
transfer matrix calculations. The observed Gonihedric model transition will 
thus be strongly influenced by crossover effects close to the endpoint of the 
standard Ising line which which is near, but not exactly at, j = —0.25. 

In such a case we would expect the correlation length to diverge as 



with v and u being the tricritical and 3-D Ising values for the correlation length 
exponent respectively, and A = j + 0.25. While CVPAM and transfer matrix 
calculations [19l[2Qj are both in agreement with this general picture the values 
obtained for the tricritical exponents and crossover exponent do not agree, the 
transfer matrix calculations finding i/ = 0.45(15) and = 0.6(2), whereas the 
CVPAM calculation finds = 1.1(1). 

In summary, while the general features of the Gonihedric model equilib- 
rium phase diagram are understood both for k ~ and k / and a plausible 
tricritical/crossover scaling scenario has been advanced to describe the ob- 
served values of the critical exponents, these have still not been determined 
with any great degree of certainty. It is also difficult to define a suitable mag- 
netic order parameter for the low-temperature phase as the flip symmetry of 
the Hamiltonian means the standard magnetisation will be zero in the ground 



(18) 



where the new feature is the critical amplitude, which is given by 



(19) 
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state. Although both low temperature expansions and CVPAM calculations 
of the free energy strongly suggest that this symmetry is broken for non-zero 
K at finite temperature, the standard magnetisation is still measured to be 
zero in the low temperature phase in Monte-Carlo simulations for k ^ Q. The 
flip symmetry appears to be unbroken for the k = (pure plaquette) case, so 
the low temperature phase for this is highly degenerate. 

A further difficulty which has been observed in attempting to extract crit- 
ical exponents via Monte Carlo simulations is the very long autocorrelation 
times which the model exhibits, as was first remarked in [5T]. These have the 
effect of smearing out the observed critical singularities and rendering the ex- 
traction of critical exponents problematic. With this in mind, we now turn to 
the non-equilibrium properties. 

3 Non-Equilibrium Behaviour, mostly by Monte Carlo 
simulations 

On heuristic grounds we would expect slow dynamics in the Gonihedric model 
since the destruction of a spin droplet will require an activation energy because 
of the weighting of any new edges in the Gonihedric Hamiltonian. Indeed, 




(a) (b) (c) 



Fig. 7. Destroying a cubic excitation, by "nibbling" away the edges has an activation 
energy in the Gonihedric model. 

when K = 1 the plaquette term drops out and one is left with 

if = -2^a,a, + i '^^'^J- (20) 

which represents a particular case of a class of models with competing nearest 
neighbour and next-nearest neighbour interactions formulated by Shore and 
Sethna specifically to investigate slow dynamics [22]. We will denote such 
models with generic coupling ratios Shore and Sethna models below. The 
Gonihedric Hamiltonian with k = 1 displays some differences in its dynamics 
compared with the generic case. 

Hand in hand with such slow dynamics and energy barriers we would also 
expect to see metastability in the model and this is backed up by both CVPAM 
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calculations and direct integration in Monte Carlo simulations to obtain the 
free energy. For k = the CVPAM calculations of the free energy show a 
region of metastability both above and below the first order transition point 
at /3 = T^^ ~ 0.55 [IS]. To calculate the free energy of the model in a Monte 
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Fig. 8. The free energy from a CVPAM calculation when k — 0. The solid line 
is the low temperature phase and the dashed dotted hne the high temperature 
paramagnetic phase. 



Carlo simulation [53] we used the following equations: 

4>IovjT 
0highT 

In the above equation C and E denote the specific heat and internal energy, 
respectively (measured using the standard formulae in the Monte Carlo sim- 
ulation), and s(oo) = ln(2) is the entropy per site at infinite temperature. 



-Tsioo)+Tj Ed{-) 



(21) 



We also find, see FiglTO] and as advertised earlier, that the lamellar low 
temperature phase has a higher free energy than the ferromagnetic phase 
when K = 1 indicating that the true low temperature phase is ferromagnetic 
in this case, although this does not appear to be realised within the timescale 
of Monte Carlo simulations carried out so far. 

Qualitative numerical experiments provide plenty of evidence for interest- 
ing dynamical behaviour in the model, particularly with respect to quenches 
[23] , In the k = case to study the evolution of a random configuration 
quenched to low temperature, one can measure the energy excess 5E{t) — 



14 



D.A. Johnston, A. Lipowski and Ranasinghe P. K. C. Malmini 
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Fig. 9. The free energy from a Monte Carlo simulation when k = 0. The transi- 
tion point and the region of metastability are both in excellent agreement with the 
CVPAM calculation of FigE] 



^ -3.6 




0.42 0.44 0.46 



0.2 0.2."^ 03 0.35 0.4 0.45 0.5 



Fig. 10. The free energy from a CVPAM calculation when k, — 1. The solid line is 
the ferromagnetic phase, the dashed dotted line the high temperature paramagnetic 
phase and the dashed line the lamellar phase 



E{t)—Eo over the ground state energy E^. One expects that the inverse of this 
quantity sets the characteristic length scale of the system l{t), which roughly 
corresponds to the average size of domains. There is convincing evidence that 
the generic behaviour in many systems with nonconservative dynamics and 
scalar order parameter (i.e., conditions which are satisfied here) l{t) increases 
asymptotically in time as l(t) ~ t" and n — 1/2. However, in some sys- 
tems l{t) is known to increase much more slowly in time, even logarithmically 
l{t) ^ log(f). These rather exceptional systems include some random (at the 
level of the Hamiltonian) systems, and also the Shore and Sethna models at 
sufficiently low temperatures. It is the energy barriers developing in these sys- 
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terns during the coarsening which cause such a slow increase of l(t), so the 
Gonihedric model would also be a good candidate for such behaviour. 

The log-log plot of l/6E{t) as a function of time for the Gonihedric model 
with K = is shown in Fig[Tl] The presented results are obtained for L = 40 
but very similar behaviour was observed for L ~ 30. From Fig llll it is clear 




1 2 log),) 3 



Fig. 11. The excess energy over the groundstate vs time during a quench to vari- 
ous temperatures. The dotted line represents generic coarsening behaviour with an 
exponent of 1/2. 

that for T — 0.75 and 1.4 the asymptotic slopes of the curves are much smaller 
than 1/2 and there is a tendency for these curves to bend downwards. Taking 
into account the absence of models with n considerably smaller than 1/2 and 
the existence of energy barriers in the Gonihedric model of basically the same 
nature as in the Shore and Sethna model, suggests that for the examined tem- 
peratures the characteristic length asymptotically increases logarithmically in 
time. It would appear that such a slow increase of l[t) takes place even for 
T = 1.65 and 1.7, but the behaviour of l{t) for these temperatures is obscured 
by the metastability effects, since before collapsing into the glassy phase the 
system remains in the liquid state for some time. 

The difference between the Gonihedric model and the Shore and Sethna 
model becomes clear when we approach the lower boundary of the metastable 
region which we roughly estimate to lie at T = Ti ^ 1.7 by increasing the tem- 
perature. In the Shore and Sethna model for temperatures below the critical 
point but above a corner-rounding transition thermal fluctuations roughen 
corners of domains and energy barriers are no longer relevant. Consequently, 
the standard coarsening dynamics with n = 1/2 is restored and the system 
rapidly evolves toward the low-temperature phase. On the contrary, in the 
Gonihedric model for 1.7 < T < 1.95 a random quench does not even evolve 
toward the low-temperature phase but remains disordered. 

The region of metastability can also be investigated through measuring 
various characteristic times by imposing different initial and boundary condi- 
tions and monitoring the evolution toward a final state. To check the decay 
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time associated with the stability of the high temperature "Uquid" phase (Tuq), 
a random initial configuration can be used and then simulated at lower tem- 
perature until the energy reaches a sufficiently small value. To estimate 
around 100 independent runs were necessary. The results for T = 1.75 are 
shown in Fig[T2] and they suggest that the escape time increases at least ex- 
ponentially with the system size, which is at first sight surprising for a model 
with finite range interactions. 

A later order of magnitude calculation [21] for the size of a critical droplet 
(one that would grow and eventually take over the whole system with the 
stable phase) in the Gonihedric model by Swift et.al. showed that the results 
discussed here were for system sizes comparable or smaller than such a droplet, 
so the system size itself was setting the scale of energy barriers, leading to the 
observed exponential growth in lifetimes. The energy barrier A to forming a 
droplet of size R would be 



A = AaR'^-^ - BSfR'' 



(22) 



in d dimensions, where a is the surface tension between stable and metastable 
phases, Sf is the bulk free energy difference between the two phases and A 
and B are determined by the geometry of the droplet. Maximising A gives 
the critical radius R* 

A{d-l)a 



R* 



(23) 



BSfd 

At low T in 3-D this gives a timescale for the nucleation of such a critical 
droplet of the form 

T - e^p{AA^a^/27B^SfT) (24) 

and inserting plausible estimates for a and Sf in the k = Gonihedric model 
gives R* = 25 for cubic droplets and r ~ 10"*^ Monte Carlo steps. F 
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Fig. 12. The size dependence of the logarithm of the escape times ruq (□) and r+_ 
(+). Calculation of ruq and t-\ — was done for T — 1.75 and T = 1.8, respectively. 



To check the stability of the low temperature "crystal" phase, one should 
measure the time (r_|_) needed for the crystal to be transformed into the liquid. 
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It would be particularly interesting to examine the size dependence of for 
1.8 < T < 1.95, i.e., for temperatures where the crystal is metastable. It was 
found, however, that this quantity increases very rapidly with the system size 
and in this temperature range it is virtually impossible to increase the system 
size beyond L — 6. The stability of this phase might be also inferred from 

measurements of the time (r-| ) needed to shrink a cubic excitation of size L. 

This technique parallels that which has already been applied to the Shore and 
Sethna model: the initial configuration has "up " spins at the boundary of the 
cube of size L + 2 (which are kept fixed) and "down " spins inside this cube. 
Simulations are performed until the magnetisation of the interior of the cube 
decays to zero and the time needed for such a run is recorded. To calculate t+_ 
at a given temperature 100 independent runs were again made. The results 

for T — 1.8 suggest that increases approximately exponentially with L, 

whereas above the corner-rounding transition in the Shore and Sethna model 

one expects L^, in which case the data in FiglT^] would bend noticeably 

downwards. 

Other analogous numerical experiments may be carried out for non-zero n 
and the tentative picture which emerges [SS] is shown in Fig[T31 The salient 

T,. 

(a) 



Tg Te 

1 I I " 

(c) 

Fig. 13. The phasehnes for the Shore-Sethna model (a), the Gonihedric model with 
a continuous transition (b) and the Gonihedric model with a first order transition 
(in particular k = 0) (c). 

features are that generic Shore Sethna models with competing nearest neigh- 
bour and next nearest neighbour interactions have logarithmically slow low 
temperature dynamics, with a transition to standard coarsening be- 
haviour above a corner rounding transition at Tcr- For the Gonihedric model 
with a continuous transition {e.g. k = I) simulations suggest that this corner 
rounding transition is pushed close to (or may even be at) the equilibrium 
critical temperature Tc. The k = Gonihedric model displays rather different 
properties. In this case there is a region of strong metastability on either side 
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of a first order equilibrium transition at T^. We tentatively identify the lower 
boundary of the region of metastability with a glass transition temperature 
Tg. We now discuss why it is tempting to do so. 

Glassy effects in simulations are usually discerned by looking at appropri- 
ate two-time quantities such as the spin-spin autocorrelation 

C{t,t^) = {^^(J^[t^)a^{t + t^)) (25) 

i 

where N = and is a waiting time before measurements are commenced. 
Both normal coarsening and glassy dynamics are expected to display strong 
waiting time dependence ( "aging" ) , but one may distinguish the type I aging 
dynamics seen in coarsening from the type II aging seen in glassy systems by 
looking at the overlap distribution function 

Q{t^+t,t^ + t) = (^^(Tr^(i + i^)af^(i + t^)) . (26) 

i 

This is measured by relaxing the system from a disordered state for a time 
tw at which point it is cloned into two sets of spins, {cr^^^} and {cr^^^}, and 
evolved with independent random numbers for a further time t. One finds 
that Q{ti^+t, tw + t) ^ a.s C{t, t^) ^ in type II aging - the heuristic idea 
being that no matter when the two systems are separated their configurations 
continue to move apart, which is suggestive of a truly rugged free energy 
landscape as found in glasses. If this were not the case, as in type I aging, one 
would find Q{tw +t,tiu +t) q for some constant q as C{t, t^) 0. For the 
K = Gonihedric Ising model type II aging is clearly seen with these criteria 
[H [551 [53] 

A further glassy feature is the behaviour of energy autocorrelation [51] 

A{t,t^)^ {E{t^)E{t + t^)) (27) 
for Tg < T < Tc- It is found that a stretched exponential form 

A{t,t^,) ^ Aoexp (- j (28) 



fits this quite well, where the relaxation time diverges as 

A 



(29) 



and the value for T* is very close to the measured value for Tg. The exponent 
P is temperature dependent and has been measure to fall in the range 0.6 — 0.8 
Further numerical cooling experiments suggest that it is possible that the 
Gonihedric model with k = may be displaying not just glassy behaviour, 
but ideal glassy behaviour ^7]. Some time ago Anderson proposed that the 
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glassy transition, which is a kinetic phenomenon, might be related in the limit 
of vanishingly small cooling rate with a thermodynamic transition [55] . It is 
found that for the Gonihedric Ising model the peak in the specific heat of 
the liquid occurs at, or very close to, the temperature at which the internal 
energy jumps under very slow cooling. The specific heat shown in Fig[T4]is a 
thermodynamic quantity. It was calculated in a (standard) quasi-equilibrium 
manner: after fixing a temperature we relaxed the system and then measured 
the variance of the internal energy. The sharp peak seen in Fig[T4] indicates a 
thermodynamic-like singularity in this model. 

On the other hand an almost vertical drop of the internal energy under 
continuous cooling, shown in Fig[T5l indicates the dynamic, glassy transition. 
The results were obtained by simulating the model under continuous cooling 

5 I ^ ^ ^ ^ 1 

4 
3 

O 

2 - 

1 ^ . 

°3 3.2 3.4 21 3-6 3.8 4 

Fig. 14. The specific heat C as a function of temperature T calculated from the 
variance of the internal energy for L — 24 (o) and L = 40 (+). At each temperature 
the system was relaxed for 10^ Monte Carlo steps and then measurements performed 
during 10'' Monte Carlo steps. 

with a constant cooling rate r and initial temperature Tq — 2.1 (Tq > Tg). 
This means that temperature as a function of time is given by T = Tg — rt. 
Calculations were performed for several system sizes L in order to ensure that 
L was sufficiently large. For example for r — 0.02 the system size L — 30 
is sufficient to obtain size-independent results but for r — 0.00002 we had to 
take L = 70. One can see that although r decreases by three decades, the zero- 
temperature energy approaches the ground-state energy very slowly and that 
lowering the cooling rate r, sharpens the transition. At first sight one might 
expect that in the limit r — > the transition becomes infinitely sharp and 
coincides exactly with thermodynamic singularity at the peak of the specific 
heat. However, as we have seen, the metastability of the liquid is a finite (but 
large) time/size effect and neither the peak nor the internal energy drop can 
be made perfectly sharp. 
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Fig. 15. The internal energy i5 as a function of temperature for (from the top) r = 
0.02, 0.002, 0.0005, 00002, 0.00005 and 0.00002. 

Anderson's idea has had a rather hmited experimental support. The main 
problem is that under slow cooling real liquids do not get trapped in the glassy 
phase but instead tend to crystallise. The reason for that is that when liquid is 
cooled below the melting point it becomes metastable and within a finite time 
due to heterogeneous or homogeneous crystal nuclcation it then crystallises. 
Only under sufficiently fast cooling can the crystal nucleation be avoided and 
the liquid trapped in the glassy state. In this context, the Gonihedric model 
with plaquette interactions appears to correspond to an almost ideal glass 
with an extremely large lifetime for the metastable state. 

A final subject which merits further investigation is the nature of the self- 
induced disorder in the glassy phase. To get some insight into the k = case it 
is possible to look the distribution of unsatisfied plaquettes in the glassy phase 
(i.e., plaquettes contributing energy above the ground state) [IS]. A random 
high-temperature sample was slowly cooled down to zero temperature and 
the final configuration was used to locate unsatisfied plaquettes and their 
spatial distribution is shown in Fig llGI For comparison we also present similar 
simulations for the k = 2 case. One can see that in both cases the energy 
from these unsatisfied plaquettes is concentrated in linear segments (around 
90% in both cases for slow cooling). Interestingly, low temperature excitations 
observed in the 3D Edwards Anderson spin glass 



where the nearest neighbour couplings {Jij} are quenched Gaussian random 
variables with zero mean and unit variance, appear to have a rather similar 
stringy structure [301 ■ It is possible that the self- induced disorder in the k = 
Gonihedric model is giving rise to similar characteristics for the excitations 
as the quenched disorder which is imposed a priori in the finite-dimensional 
spin glass. 




(30) 
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Fig. 16. A distribution of unsatisfied plaquettes in the zero-temperature glassy 
phase. Simulations were made for the system size L = 100. 



A further point of contact with spin glasses is that a spin glass version of 
the K = Hamiltonian (still on a 2>D cubic lattice) 



(31) 



where the Jn are now quenched random variables, was also found to dis- 
play numerous common properties with structural glasses '31'. Many of the 
observed numerical features seen in simulations of this model, such as the 
stretched exponential correlations and divergence of the autocorrelation time 
are extremely similar to those discussed here for the k = Gonihedric model. 
The inference from this might be that the additional quenched disorder aris- 
ing from the Jn may be gratuitous, since the plaquette Hamiltonian on its 
own is enough to give the same, or very similar, behaviour. 

Self-induced disorder has been posited before as a mechanism for glassy be- 
haviour in various deterministic mean-field models. One well-known example 
is the Hamiltonian originally formulated by Golay [32] 



H 



7^2 

•Jq \ ' 

O.N 2^ 



2N 



k=l 



N-k 

E 



(32) 



in the context of obtaining binary sequences with small autocorrelations. It 
was shown that many of the properties of the deterministic Hamiltonian in 
eau.(|32p were reproduced by the disordered Hamiltonian 



H 



1 

'2N 



N-l 

E 

fc=i 



N 

E^, 

2J = 1 



(fc). 



(33) 



where the J^^'^ were random connectivity matrices |33j . The contrasting fea- 
ture for the Gonihedric model is that the self-induced disorder is appearing 
with finite range couplings. 

The non-equilibrium, dynamical properties of the Gonihedric Ising model 
thus display many interesting features both for k = and k ^ 0. For the 
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case of K = (the purely plaquette action) there are clear parallels with the 
observed behaviour of real glass formers, so it offers an appealing lattice model 
for investigations of glassy dynamics by virtue of its simplicity. Indeed, there 
may also be similarities with finite dimensional spin glasses, as we have seen 
above. 



4 Variations on the Glassy and Gonihedric Themes 

There are a plethora of other lattice models which can be formulated to display 
glassy properties [MlISS]. In many cases, these have trivial static Hamiltonians 
of the form 

N 

H = Y,n,, (34) 

i=l 

where the = 0, 1 are defined on each site of a cubic lattice of linear size L, 
which has periodic boundary conditions and can be though of as labelling the 
(im)mobility of a site. The interesting (glassy) behaviour is induced in these 
cases via constrained dynamics. 

Other models may possess a non-trivial Hamiltonian, but a trivial equi- 
librium phase diagram. Nonetheless, the dynamical properties of such models 
may still be of interest. The 2D Gonihedric Ising model is one such case 



^gon^/^ ^^^^ ^^(^3 + '^ ^^ 3 ^ ^ ^r^^J 



1 - K 



where the form is similar to the 31? model, but the relative weights of the 
terms differs. This is constructed to weight the corners of 2D spin clusters, 
rather than their edge length, which is the effect of the standard 2D Ising 
action. Via a mapping to the 8-vertex model it is possible to show that the 
2D Gonihedric Hamiltonian does not display a static phase transition. Just 
as for the trivial Hamiltonian models with constraints, the dynamics displays 
many interesting features [36j . 

We have concentrated exclusively on Ising spins here. It is also possible to 
define a 3D Gonihedric Potts models 

H = -2^Y. (2^-.,-.+. - 1) + f E (2'5...-.+.+. - 1) 

a;,/27^i^,/2>0,i/>0 

(2(5c,-x,(Tx + ^'^crx+,, ,crx+j,+i, ~ 1) ' 

(2'5ax,(Tx+^+.'5^x+.,(Tx+^ - 1) (35) 

where the spin variables {a} now take on q values [37]. In this case the model 
appears to retain a first order equilibrium transition for larger k than with 
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Ising spins, but the dynamical properties of the model remain to be investi- 
gated. 

The idea that multi-spin, specifically plaquette, interactions give rise to 
self-induced disorder might also be applied to models with continuous spins 
such as gauge glasses. The standard 3D gauge glass Hamiltonian is of the 
form [38] 

H = - ^ cos (61, - Oj + A,j) (36) 

<i,3> 

where the Aij are quenched, random variables uniformly distributed in [0, 27r]. 
A plaquette interaction of the form 

H = - J2 cosie,^e, + 9k-ei) (37) 

lt,j,kd] 

would presumably give rise to a similar sort of frustration as the quenched 
disordered {Aij}. 

A final avenue for further research might be the dual model to the k — 
Gonihedric action, which may be written in several equivalent forms including 

Hduai = ^''(^^'^(e +x)+ ^"(c)A"(e AHmH^ + c) (38) 

where = (1, 1, -1, -1), A'l = (1, -1, 1, -1) and A'^ = (1, -1, -1, 1) are one 
dimensional irreducible representations of the fourth order Abelian group and 
^, 1], ^ are unit vectors in the dual lattice. The spins may also be considered as 
Ising (±1) spins if we set A'^ = A^ A^, which gives the following Hamiltonian 

Hduai = E (^)^'' + ^) + (^)^'' + + (^)^' (^)^'' + c)^" (e + c) ■ 

(39) 

This is recognizable as a strongly anisotropically coupled Ashkin- Teller model. 
It is possible that consideration of this dual model might shed some light on 
the definition of an order parameter in the original Hamiltonian and also 
suggests an intriguing link with anisotropic scaling properties. 

5 Endpiece and Acknowledgements 

Gonihedric 31? Ising models have had an interesting past and display novel 
properties both in equilbrium and dynamically. The k = (plaquette) Go- 
nihedric Hamiltonian appears to provide a simple lattice model with similar 
behaviour to structural glass formers, which also has some connections with 
spin glasses in spite of the absence of quenched disorder. Even after some fairly 
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extensive investigation by Monte Carlo simulations and other means over re- 
cent years, many features of the models still remain in need of clarification. 
They may well have an interesting future ahead of them too. 
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